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Abstract 

The bouncing cosmology provides a successful solution of the cosmological singu- 
larity problem. In this paper, we study the bouncing behavior of a single scalar field 
model with tachyon field non-minimally coupled to itself, its derivative and to the 
curvature. By utilizing the numerical calculations we will show that the bouncing so- 
lution can appear in the universe dominated by such a quintom matter with equation 
of state crossing the phantom divide line. We also investigate the classical stability of 
our model using the phase velocity of the homogeneous perturbations of the tachyon 
scalar field. 
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1 Introduction 



Cosmological observations have revealed that our universe has experienced a phase of ac- 
celerated expansion in the early-time known as inflation [1]. Although inflation has great 
successes such as solving the flatness, horizon and homogeneity problems [2], but it suffers 
from the initial singularity problem and hence cannot give complete descriptions of the early 
universe. 

A successful solution of the cosmological singularity problem may be provided by non- 
singular bouncing cosmology which has been proposed a long time ago [3]. In order to 
circumvent the Big-Bang singularity problem one needs to choose one of the following ap- 
proaches: 1) modified gravity [4], 2) using a kind of matter field that violates the null energy 
condition (see [5] for review). Through the modified gravity approach various scenarios have 
been constructed such as the brane world models [6], Pre-Big-Bang [7] and the Ekpyrotic 
models [8], f(T) gravity in which gravity is described by an arbitrary function of the torsion 
scalar [9], higher derivative gravity actions [10], non-relativistic gravitational actions [11, 12] 
or loop quantum cosmology [13]. Bouncing cosmologies can also be inspired by string theory 
when the string gas cosmology [14, 15] embed in a bouncing universe [16]. 
Furthermore, Bouncing cosmology for a homogeneous and isotropic universe can lead to the 
matter bounce scenario. Non-conventional fluids [17, 18], non-minimally coupled scalar fields 
to gravity [19] and ghost condensates [20-22] are belong to this sort of bouncing cosmology. 
The simplest possibility to obtain a bounce is a massive scalar field in the presence of a pos- 
itive scalar curvature [23] and the other type of the matter bounce cosmology is the quintom 
bounce model [24, 25]. 

Moreover, various observational data such as SNe la Gold data set [26] confirmed that the ef- 
fective equation of state (EoS) parameter u (the ratio of the effective pressure of the universe 
to the effective energy density of it) cross —1, namely, the phantom divide line, currently or 
in the past. 

The scenarios with phantom behavior are designed to understand the nature of dark energy 
with bj across —1. To realize a viable scenario with phantom behavior one needs to introduce 
extra degree of freedom to the conventional theory with a single fluid or a single scalar field. 
The fact that the phantom behavior can be achieved without ghosts in scalar-tensor gravity 
was first shown in [27]. Also, quintom scenario which realizes the phantom behavior first 
introduced in Ref. [28] with two scalar fields (quintessence and phantom). Another attempts 
for constructing a model that shows phantom behavior are as follows: scalar field model with 
non-linear kinetic terms [29] or a non-linear higher- derivative one [30], braneworld models 
[31], phantom coupled to dark matter with an appropriate coupling [32], string inspired 
models [33], non-local gravity [34], modified gravity models [35] and also non-minimally cou- 
pled scalar field models in which scalar field couples with scalar curvature, Gauss-Bonnet 
invariant or modified f(R) gravity [36-38] (for a detailed review, see [39]). Crossing of the 
phantom divide can also be realized with single imperfect fluid [40] or by a constrained single 
degree of freedom dust like fluids [41]. 

Moreover, non-minimal couplings are generated by quantum corrections to the scalar field 
theory and they are essential for the renormalizability of the scalar field theory in curved 
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space (see [42] and references therein). One can extend the non-minimally coupled scalar 
tensor theories, allowing for non-minimal coupling between the derivatives of the scalar fields 
and the curvature [43]. A model with non-minimal derivative coupling was proposed in [43- 
45] and interesting cosmological behaviors of such a model in inflationary cosmology [46], 
quintessence and phantom cosmology [47, 48], asymptotic solutions and restrictions on the 
coupling parameter [49] have been widely studied in the literature. General non-minimal 
coupling of a scalar field and kinetic term to the curvature as a source of dark energy has 
been analyzed in [50]. Also, non- minimal coupling of modified gravity and kinetic part of 
Lagrangian of a massless scalar field has been investigated in [51]. It has been shown that 
inflation and late-time cosmic acceleration of the universe can be realized in such a model. 
In this paper we consider a model with an explicit coupling between the scalar field, the 
derivative of the scalar field and the curvature and study the bouncing solution in such a 
model. We are interested in our analysis to the case of tachyon scalar field. The tachyon 
field in the world volume theory of the open string stretched between a D-brane and an anti- 
D-brane or a non-BPS D-brane plays the role of scalar field in the context of string theory 
[52]. What distinguishes the tachyon Lagrangian from the standard Klein-Gordan form for 
scalar field is that the tachyon action has a non-standard type namely, Dirac-Born-Infeld 
form [53]. Our motivation for investigating a model with non-minimal derivative coupling 
and tachyon scalar field is coming from a fundamental theory such as string/ superstring the- 
ory and it may provide a possible approach to quantum gravity from a perturbative point 
of view [54-56] . 

An outline of the present work is as follows: In section 2 we introduce our model in which 
the tachyon field plays the role of scalar field and the non-minimal coupling between scalar 
field, the derivative of scalar field and Einstein tensor is also present in the action. Then, we 
derive the field equations as well as the energy density and pressure of the model in order to 
study the bouncing behavior of our model. We obtain the condition required for bouncing 
solution of our model and then we show that such a condition can be satisfied numerically. 
In Section 3 we discuss the classical stability of our model. Section 4 is devoted to our 
conclusions. 

2 Bouncing behavior of non-minimal derivative cou- 
pling of tachyon gravity 

In this section we review the field equations of tachyon gravity with non-minimal derivative 
coupling and provide the cosmological equations in a Friedmann-Robertson- Walker (FRW) 
universe. Then, we obtain the necessary conditions required for a successful bounce in such 
a model. Our starting point is the following Born-Infeld type action for tachyon field with 
non-minimal derivative coupling and also with itself, 
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where k 2 = 8tcG = y4- while G is a bare gravitational constant and Mpi is a reduced Planck 

mass, A is the cosmological constant and V(4>) is the tachyon potential which is bounded 
and reaching its minimum asymptotically. f{4>) is a general function of the tachyon field 
4> and £ is coupling constant. The models of kind (1) with non-minimal coupling between 
derivatives of a scalar field and curvature are the extension of scalar-tensor theories. Such 
a non- minimal coupling may appear in some Kaluza-Klein theories [57, 58]. In Ref. [43], 
Amendola has considered a model with non-minimal coupling between derivative of scalar 
field and the Ricci scalar, ^Rd^d^cf), and by using generalized slow-roll approximations, he 
has obtained some inflationary solutions of this model. 

A general model containing two derivative coupling terms ^iRd^d^cp and t^R^d* 1 §d v <fi, has 
been studied in [44, 45]. It was shown in [47] that field equations of this theory are of third 
order in g^ v and (f), but in the special case where — 2£i = £2 = £ the order of equations are 
reduced to the second order. This particular choice of £1 and £2 leads to the non-minimal 
coupling between derivative of scalar field and the Einstein tensor, ^G tiu d ll (pd u (j). Sushkov in 
[47] has obtained the exact cosmological solutions of this theory and he has concluded that 
such a model is able to explain a quasi-de sitter phase as well as an exit from it without any 
fine-tuned potential. 

Let us now present the cosmological equations. Variation of action (1) with respect to the 
metric tensor g^ u gives, 

R^ - \g^R + = k 2 (T^ + £t;J , (2) 

where 

T„ = VMp£!4- Slm u), (3) 

and 

T;, = R(V,<PV U <P) - 4V 7 0V (/ ^ ) + G^(V0) 2 - 2i^ A WV A <£ - 2V M WV„V 7 

+ 2 V M V v (j>n<j) + 9tlv ( V V A 0V 7 V A -(□</>) 2 + 2W X V 7 0V A 0) , (4) 
where u = yT + V^V^. 

One can obtain the scalar field equation of motion by varying (1) with respect to 0, 

V„(^P) - d -m u - m)G ^^ + ^G^m = 0. (5) 
For a flat FRW metric, 

ds 2 = -dt 2 + a 2 (t)(dr 2 + rW), (6) 
the components of the Ricci tensor R^ v and the Ricci scalar R are given by 

R 00 = -3(H + H 2 ), Rij = a 2 (t)(H + 3H 2 )5 ij , R = Q(H + 2H 2 ) , (7) 
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where H = is the Hubble parameter and a(t) is the scale factor. The scalar field equation 
of motion for a homogeneous (f> in FRW background (6) takes the following form 

6 • 1 dV 

— ?-t- + 3H(f) + 



V(<j)) d(j> 



+ V ^ (3£tf 2 (2/(0)0 + ^0 2 ) + 18£tf 3 /(</>)0 + 12tHHf(<f>)jA = 0. (8) 
From equation (2) the energy density and pressure are as follows respectively, 

p = 4^= + 9^ 2 /(0)0 a + 4. ( 9 ) 



and 



P = -Vffly/l- <P 2 - £{3H 2 + 2iT)/(0)0 2 - 2^(2/(0)# + |^ 3 ) - A (i ) 



Friedmann equation is also as follows, 



2 _ K 2 ( V{(j>) , ^^^^^ , A 



h 2 = y[-j^= + H&MW) + 3 • (ii) 



Moreover, a bouncing universe starts with an initial contraction to a non-vanishing minimal 
radius, then subsequent an expanding phase [24]. During the contracting phase, the scale 
factor a(t) is decreasing, i.e., a < 0, and in the expanding phase we have a > 0. At the 
bouncing point, a = 0, and around this point d > for a period of time. Equivalently in the 
bouncing cosmology the Hubble parameter H runs across zero from H < to H > and 
H = at the bouncing point. A successful bounce requires around this point, 



K 2 , _ 



= -—{p + P) = --p{l +u)>0. (12) 

As it is clear from equation (12), a successful non-singular bouncing cosmology violates the 
null energy condition |(NEC) around the bouncing point. In addition, in order to achieve 
the hot Big-Bang era after the bouncing, the EoS of the universe, u, must cross the so called 
phantom divide line u — — 1. 

A salient feature of our model (1) is that it realizes the phantom divide line crossing [59]. 
In Ref. [59] we have presented the conditions required for such a crossing and numerically 
shown that this model realizes crossing of the u = — 1 for a special case of the Hubble 
parameter H = with constant h . In the present work we show that crossing over —1 can 
take place in our model for a general Hubble parameter obtained from equation (11). 



3 The NEC for a null vector k^, implies T^k^k v > which for a perfect fluid reads p + P > 0. 
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To study the bouncing behavior of the present model we use the equation (12). From 
equations (9), (10) and (12) one can see that a successful bounce requires: 



' ^H 2 f{4>)4> 2 - 2iHf{4>)4> 2 - 2£H (2/(0)# + |^ 3 ) < 0. (13) 



1 - 



We will show below that our model can easily fulfilled the above bouncing condition. Note 
that in the numerical study on the bouncing solution, we have used the following state for 
Hubble parameter obtained from equation (11): 

In figures 1 and 2, we show the bouncing solution for two different potentials. Figure 1 shows 
such a bouncing solution for V(4>) = Voe~°"^ with constant a. One can see from this figure 
that the model predicts crossing of the uj — — 1, which gives rise to a possible inflationary 
phase after the bounce. Also, the Hubble parameter evolution as a function of cosmic time 
has been shown in figure 1. One can read from this figure that the Hubble parameter H 
running across zero at t — which have choosed it, as the bouncing point. In figure 2, 
we consider the model with potential V((f>) = ^ and then show another example of the 
bouncing behavior of the model. This figure show the crossing over —1 for the EoS again. 
In the numerical calculations we have used the function = bcf> n with constants b and n. 




Figure 1 : Plots of the evolution of the EoS and the Hubble parameter versus t for the potential 
V{4>) = Vbe - "^ , <fi = fat, f(4>) = b<j) n and negative cosmological constant A = —0.05, (with £ = 10, 
b = 1, n = 8, V = 0.04, 0o = 0.6 and a = 5). 
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Figure 2: Plots of the evolution of the EoS and the Hubble parameter versus t for the potential 
V{4>) = % <j) = 4> t, f(4>) = b4> n and A = -0.05, (with £ = 10, b = 1, n = 8, V = 0.04 and 
00 = 0.6). 

3 Stability of the model 

Now we discuss on the stability of the model. The sound speed, C 2 , is the function appearing 
before spatial gradients in the scalar field equation of motion and for perturbations around 
homogeneous solutions, it is a function of time. The sound speed express the phase veloc- 
ity of the inhomogeneous perturbations of the tachyon field [60]. The absence of gradient 
instabilities requires C 2 > 0, where 

Cl = j, (15) 

where a prim denotes derivative with respect to 2 . By using equations (9), (10) and (14), 
the sound speed parameter reads, 



3(J,-27£/(0)0 2 + 18£/(0)) 



3 ^ W(l -0 2 )# 



9 , 4 4f , 70 2 2ie/(0)0V fi . fU uri ^ 3, 

"w + 7^ 771 + 77777 177 ~ "7T77 177" + 6 £/W 1 ~ & ) ~ o< 

fc 2 4 d0 4P(1 - 2 ) 4(1 - 2 ) 2 



135£ 2 / 2 (0)0 3 L^2 n l2^df M 4> ~tfU\jJk 

— ( 3£0 (1 - )— + - - 3£/(0)00 



2(^-9e/(0)0 2 ) 5 (i 



+6£/(0)0(l - 2 ) ] + 27^ 2 / 2 (0)0 2 (l - 2 ) + 3^/(0) (1 - 9£/(0)0 2 ) (1 - ^0 2 ) 



7 



2\/(l-</> 2 ) 



9e/(0)^ 2 ) 3 (i 



«i(^| + /(^ + 2(l-^) 



-9? 2 /(0)<% 



9e/(^ 2 )Mi-0 2 ) ! (^ + 2 



+ 



4 



27e/(0)0 2 (l - fi) 
2(J,-9£/(0)0 2 ) 



(16) 



This relation is very complicated to find exact analytical conditions for the stability of 
the model. So, we perform the numerical computations to show the evolution of the C 2 . In 
figure 3, we have plotted the C 2 for the models considered in this paper for the numerical 
calculations shown in figures 1 and 2. From this figure, we can see that the sound speed 
parameter is positive throughout the bouncing phase. 
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Figure 3: Plots of the sound speeds versus t, (left for the potential V(4>) = V$e a< ^ and right for 
the potential V{4>) = ^) , $ = fat and f((f>) = b<j> n , (with f = 10, b = 1, n = 8, V = 4, fa = 0.5 
and a = 5). 



Moreover, we should expect that the coefficient of in equation (8), is responsible for 
the presence or absence of ghosts. We show this coefficient by D (for more details see [61]). 
In order to ensure that the perturbation of </> are not ghosts, one needs D > 0, where D is 
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as follows, 



D 



V(<j>) 



+ 6^ 2 /(0)- 



(17) 



(1-02): 

The evolution of D have been shown in figure 4 for two different potentials. It is clear from 
this figure that the function D is positive during the bouncing phase. 



D 




D 




Figure 4: Plots of the D versus t, (left for the potential V{4>) = V$e and right for the potential 
V(4>) = , (/> = fat and /(0) = b(f) n , (with £ = 10, b = 1, n = 8, Vq = 4, (f) = 0.5 and a = 5). 

Now we study the tensor perturbations of our model. We use the procedure developed in 
[62] for most general scalar-tensor theories with second-order field equations. It has been 
shown that the quadratic action for the tensor perturbations is as, 



S T = - I dt d 3 x a 3 



(18) 



are 



hij is a tensor perturbation satisfying = = hijj and in case of our model Gt and Ft 

G -<i-iWf), (10) 



(20) 



In order to avoid ghost and gradient instabilities in the tensor sector we require the conditions 
F T > and G T > 0. 
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We find the power spectrum of the primordial tensor perturbation as follows, 



Pt = 87t^^ 



-ky T =l 



(21) 



-\ TT-dt and k is the Fourier wavenumber. Note that the sound horizon 



where dyr 

crossing occurs when k 2 ~ \. Also in (21) we have 



here v T = 2 _ 3 ~l^l gr , where 



V{v T ) 



r(i; 



: (l-e-- + -)- 

v 2 2 



(22) 



9t 



h 



«2W + 



(23) 



(24) 



and e = —jp that takes the following form for our model 



2 (^ + £)V- 9 ^^ 2 ) 5 

Finally, the tensor spectral index is given by 

tit = 3 — 2vt- 



(25) 



(26) 



Thus, one can obtain a blue or red spectrum of tensor perturbations by choosing different 
values of coupling parameter £ in tachyonic non-minimal derivative coupling model. 
Let us focus on scalar perturbations by putting hij = 0. The quadratic action for curvature 
perturbations is as follows, 



S, 



dt d 3 x a 3 



GsC 2 



■(vc) 5 



(27) 



here C is scalar perturbation in perturbed metric (for more details see [62]) also Gs and F$ 
are given by, 

G s = ^G 2 r + 3G T , (28) 
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where H and G in the case of our model are, 



2(1 



3H' 



6^ 



(29) 



(30) 



(31) 



By substituting above equations as well as (19) and (20) into Gs and F$, one can obtain, 



G s 



i + 3^ 2 «) \2^(i - + s^«) 



TO)\ 



3 J 



i + 3^2^ 



# 7 



# 2 ^ + 3^0 2 



+ 



-20£ 



d f(<t>)\(^df((t>) 



-) 



20- 



12^7(0) \ fl ,y 



, (32) 



W)\ 2 



• (33) 



As it is mentioned in [62] , the analysis of the curvature perturbation hereafter is completely 
parallel to that of the tensor perturbation. So, ghost and gradient instabilities are avoided 
as long as Fs > and Gs > 0. 

The power spectrum of the primordial curvature perturbation is given by, 



Pr 



IsGl H 2 



2 f Hit 2 



-ky s =l 



(34) 



where dys = -\/fr-dt and k is the Fourier wavenumber. Note that the sound horizon 



a V G S 

crossing occurs when k 2 ~ \. Also in (34) we have 



Is 



n 



fs_ £te\ 
2 + 2 h 



(35) 



here we have defined vs 



3 ~ € +9s wri ere 

2-2e-/ s + 9s ' Wnt ^ 



9s 



fs 



G< 



HG S 

F s 
' HFs 
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(36) 
(37) 



The spectral index is, 



n s - 1 



3-2^ 



(38) 



The tensor-to-scalar ratio is given by, 



r = 




) 



(39) 



4 Conclusion 

In this work we investigated the bouncing solution in the universe dominated by the quintom 
matter. We analyzed the possibility of obtaining a cosmological bounce in a model which 
contains a scalar field with non-minimal derivative coupling to Einstein tensor and itself 
too. We used the numerical methods to study the bouncing behavior in this setup where 
the tachyon field played the role of scalar field. We obtained the bouncing condition as 
equation (13), then by considering a couple example for potential of the tachyon scalar field 
we showed that the above mentioned condition can be satisfied. After that we considered 
the conditions required for the classical stability of our model. These conditions obtained 
as equations (16) and (17) which reflect the positivity of the sound speed and D. Our 
numerical results for such a classical stability plotted in figures 3 and 4 respectively. Finally, 
we investigated the tensor perturbations of our model using the method of Ref. [62]. 
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